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ABSTRACT 

Pinning of superfluid vortices to the nuclear lattice of the inner crust of a neutron 
star supports a velocity difference between the superfluid and the solid as the star 
spins down. Under the Magnus force that arises on the vortex lattice, vortices un- 
dergo vortex creep through thermal activation or quantum tunneling. We examine the 
hydrodynamic stability of this situation. Vortex creep introduces two low-frequency 
modes, one of which is unstable above a critical wavenumber for any non-zero flow 
velocity of the superfluid with respect to the solid. For typical pinning parameters 
of the inner crust, the superfluid flow is unstable over length scales ^ 10 m and over 
timescales as fast as months. The vortex lattice could degenerate into a tangle, and the 
superfluid flow could become turbulent. Unexpectedly large dissipation would suppress 
this instability. 
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1 INTRODUCTION 

The dynamics of the superfluid interior of a neutron star is central to understanding a variety of phenomena that includes 
observed spin glitches, stochastic spin variations and thermal evolution, as well as possible precession and r-modes. In this 
connecti on, the possib l e importan ce of hydrodynamic instabilities in neutron stars has become a question of considerable 
interest. 



m, the possib l e imp ortance ot hydrodynamic instabilities m neutron stars has become a question ol considerable 
Peralta et all |2005l , bood ) have shown that differential rotation in the core, resulting from a spin glitch or possibly 



causing it, drives an Ekma n flow along the rotat ion axis than can excite a variant of the "Glaberson-Donnelly" counterflow 
instability in liquid helium (jGlaberson et al . 1974); transitions between laminar flow and fully-developed turbulence could drive 



spin glitches. This instability could also be excited in precessing neutron stars ( Glampedakis et al. 20081: van Hoven fc Levin 



2008). Unstable shear layers jPeralta fc Melatoslliouil) and r-mode instabilities (jGlampedakis fc Anderssonll2009l ) in the outer 
core may also play a role in glitches. 

From the standpoint of building a realistic theory of neutron star seismology with which to interpret observations, it is 
important to identify hydrodynamic instabilities of possible relevance. The possibility of turbulent instabilities in the neutron 
star inner crust, the region from the neutron drip density to about half nuclear saturation density, has received little attention 
in this regard. Here the vortices that th read the rotating superfluid are predicted to interact with nuclei with energies of ~ 1 — 5 
MeV per nucleus in the denser regions ( Alpar 1977 ; Epstein fc Bavm 19881 : Donati fc Pizzocherol 20061 : Avogadro et al. 2007 ) . 
Recent wor k has shown that this interaction will pin vortices to nuclei, regardless of the details of the pinning potential 
( Linkll2009h . As the star spins down, the differential velocity between the superfluid and the pinned vortices approaches 
the crit ical value at which the hydrodynamic lift force on vortices, the Magnus force, would unpin them. As suggested long 
ago by lAnderson fc Itohl (|l975i ). the spin glitches seen in neutron stars could arise from large-scale vortex unpinning from 
nuclei, whe rein the threshold for pinning is ex ceeded. Below the critica l velocity, pinned vortices slowly creep through thermal 
activation 1 Alpar et al. 1984 : Link et al. 19931 ) or quantum tunneling ( Link et al. 19931 ). driven by the Magnus force. Here we 
demonstrate the existence of a hydrodynamic instability related to the vortex creep process that could grow over timescales 
as short as months. 
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In the next section, we describe vortex pinning in the inner crust. In §3 we give the stability analysis. In §4, we discuss 
hydrodynamic wave solutions in the case of no background flow. In §5, we describe the hydrodynamic instability that arises 
when vortices move slowly through the nuclear lattice. In §6, we calculate the vortex mobility, which we apply to obtain 
the growth rate of the instability. We conclude with a discussion of the possibility that the inner-crust superfluid becomes 
turbulent. 



2 VORTEX PINNING 

First we calculate the critical velocity between the superfluid and the crust that can be sustained by vortex pinning. We will 
use these results in §6 to calculate the vortex mobility and the growth rate of the instability. 

Vortex pinning fixes the local superfluid velocity in the laboratory frame. As the crust spins down, a velocity difference 
v between the pinned vortices and the superfluid develops. The Magnus force per unit length of vortex is 



fn 



pKV, 



(1) 



where p is the superfluid mass density, k = h/m is the quantum of vorticity, and m is twice the neutron mass. Let F p be the 
characteristic force of the vortex- nucleus potential. Above a critical velocity difference v c , the Magnus force will exceed the 
pinning force, and vortex pinning is not possible. If a vortex could bend to intersect nuclei of average spacing a, the critical 
velocity difference v c would be given by 



pKVcCl = F p . 



(2) 



A vortex has a large self energy (tension) that typically prevents it from bending over a length scale a. If the tension 
were infinite, a vortex could not pin at all, since the vortex would remain straight and the forces from nuclei that surround 
the vortex would cancel on average. For finite tension, the vortex can bend over a length l p > a, and the critical velocity is 
given instead by 



pKvJp = Fp, 

giving a critical velocity 



Fp_ 

pKa 



(3) 



(4) 



Tension lowers the critical velocity by a factor a/l p . To calculate l p , let r v (z) be a vector in the x — y plane that gives the 
shape of the pinned vortex. The energy of a static vortex in a pinning field V(r v ), in the absence of an ambient superfluid 
flow, is 



Ev — 



dz 



dr v (z) 



dz 



+ V{r v ) 



(5) 



where T v is vortex tension, typically I MeV fm 1 . On average, over a length l p the vortex bends by an amount 5r v to intersect 
one nucleus in a volume lpir(8r v ) 2 . The quantities l p and 5r v are therefore related by 



a 3 l p n(5r v ) 2 = 1. 



The energy of the vortex per unit length, from eq. 1(5). is approximately 



Ev 



2 In 



(6) 



(7) 



p 'p ip 

where E p is the interaction energy between a vortex and a single nucleus, typically ~ I MeV. Contributions to the potential 
by nuclei that the vortex does not intersect have been ignored; these contributions will largely cancel. Minimization of E v /l p 
with respect to l p , using eq. l[6). gives 

1/2 

(8) 



3aT v 



2tcE, 



The vortex ten sion T v is due mainly to the kinetic energy per unit length of vortex due to circulation abut the vortex, and 
takes the form |Thomsonlll88ol ; iFetterll 19671) . 



T v = (0.116- In k v £), 

47T 

where £ is the radius of the vortex core and k v is the characteristic bending wavenumber, k v = n/2l p . 

For typical conditions of the inner crust, the ratio l p /a is much larger than unity. At a density p — 5 x 10 1 " g cm 



(9) 



the 
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lattice spacing is a ~ 50 fm and the radius of the vortex core is £ ~ 10 fm. For E p = 1 MeV, simultaneous solution of eqs. 
(O and ([9| gives l p ~ 9a. The ratio l p /a increases for weaker pinning. For example, for E v = 0.1 MeV, the pinning length 
becomes l p ~ 32a. For E p = 10 MeV, unrealistically large according to recent calculations, l p — 2a. 
Combining eq. (|SJ) with eq. @ gives the critical velocity, modified by vortex tension, 



F p ( a\ E„ / 2E„\ 1/2 



p , (10) 
pKa Up I p/ca£ \ 3aT v J ' 

where we have tak en F v = E v /(,. Eq. (|10|l was found in the numerical simulations of the dynamics of an isolated vortex in a 
random potential 1 Link! 2009 ). This equation shows that pinning is weakened by vortex tension. For E p = 1 MeV and £ = 10 



fm, the critical velocity is v c — 4 x 10 5 cm s _1 .0 

The corresponding differential angular velocity between the superfluid and the crust is as large as ~ 1 rad s _1 , but still 
much less than the angular velocity of the star when the superfluid condensed. The relative flow between the superfluid and 
the crust will thus be close to or comparable the local critical velocity in regions where there is pinning. We now examine the 
stability of this differentially-rotating state. 



3 PERTURBATION ANALYSIS 

The problem of th e coup led dynamics of the superfluid and vortex lattice can be studied using the hydrodynamic theory of 
Bavm fc Chandler! l|l983h which accounts for vortex degrees of freedom. The local quantities of fluid velocity, vortex density, 



and vortex velocity are averaged over a length scale that is large compared to the inter- vortex spacing l v ; the theory is valid 
for wavenumbers that satisfy kl v « 1. We treat the superfluid as a single-component fluid at zero temperature, and ignore 
dissipation in the bulk fluid and the small effects of vortex inertia. These approximations are justified in a typical neutron 
star, for which the temperature of the inner crust is much less than the condensation temperature of the superfluid. We also 
treat the crust as infinitely rigid and ignore local shear deformations, an approximation that will be justified below. The 
motion of the superfluid does not couple to the electrons, so electron viscosity is not relevant. Magnetic fields are not relevant 
either, as they do not interact with the vortices of the inner crust. 

We will consider only shear modes in the superfluid, so that the flow velocity v(r,t) is divergence-free. The rotation axis 
lies along z, and r v (r,t) denotes the continuum vortex displacement vec tor, with components in the x — y plane only. The 
equations of motion in the laboratory frame are (|Bavm fc Chandlerll 1983^ 1 



V-u = (11) 
Ov 

— +v Vv = -Vfi-Vcp- a et /p + f/p (12) 

poj x (v - ^) = -<r e i(r v ) + f. (13) 

where uj = V x v is the vorticity due to the existence of vortices in the fluid, p is the chemical potential, (f> is the gravitational 
potential, a e i/p is the elastic force per unit volume that arises from bending of the vortex lattice, and f/p is the force per 
unit volume exerted on the fluid by the normal matter. The elastic force is 

cT.i/p = -4 [2V_l(V ■ r v ) - Vlr„] + c v -^- (14) 

where denotes a der ivative with components in the x — y plane only. Here ct = (fif2/4m) 1//2 is the Tkachenko wave 
speed |Tkachenkolll966al lr}|. and Q, is the spin rate of the superfluid. The quantity Cy = (hQ./2rn) ln(Q. c /Q) is related to wave 



propagation along the rotation axis; Q c = h/{\f^m^ 2 ) for a triangular vortex lattice. For a typical neutron star rotation rate 
of SI = 100 rad s" 2 , c T = 0.09 cm s 1 and cv — 9ct- The areal density of vortices in the x — y plane is l v 2 = 2mQ./h for a 
uniform vortex lattice; hence, the requirement that kl v « 1 is equivalent to kcr « fl. 

Eq. (|13|l is an expression of balance of the Magnus force, the elastic force of the deformed vortex lattice, and the force 
exerted on the fluid by the normal matter. If the vortex array is perfectly pinned to the normal matter of the inner crust 
moving at velocity v n , so that dr v /dt = v n , the force is 

/ = /)ii)X(»-ti n ) + (rei(r,). (15) 

For imperfect pinning, the Magnus force and elastic force drive vortex motion with respect to the normal matter. For imperfect 
pinning, the force above can be generalized as 

1 A value of v c as large as ~ 10 7 cm s^ 1 was estimated in assuming E p = 5 MeV at p = 10 g cm , for £ ~ a ~ 70 fm. 

This number is a generous upper limit. 
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f = P'puj x (v - v n ) + Pp& x (w x {v - v n }) + (1 - 7) er et (r v ). (16) 
The first two terms of this force are present in the mutual friction force introduced bv lHall fc Vinen jl956h . We empha 



rasize 

the generality of the force law of eq. I|16p . The first two terms represent the force exerted on the fluid by vortices that are 
moving with respect to the normal matter; the first term corresponds to the force transverse to the vortex motion, while the 
second term corresponds to the force parallel to the vortex motion. The coefficients a and /3 can be calculated using a specific 
theory of vortex mobility. The third term accounts for the the contribution to the force that arises from local vortex bending. 
If the vortex lattice is locally undeformed (cr e i = 0), the vortex velocity from eqs. (|16jl and (|1 3 p is 

= v n + a(v - v n ) - /3uj x (v - v n ), (17) 

where a = 1 — /?'. Imperfect pinning, that is, "vortex creep", corresponds to a « 1 and /3 « 1. We refer to a, /3, and 7 as 
the "pinning coefficients". Perfect pinning corresponds to the limit a = /3 = 7 = 0, while no pinning (/ = 0) corresponds to 
a = 7 = 1 and ,5 = 0. Vortices move with a component along v — v„, so that < a ^ 1. The energy dissipation rate per unit 
volume is determined by /3, which must be positive to give local entropy production. 

Vortex creep could be a low-drag process, with /3 << a, or a high-drag process, with /3 >> a. In much previous work on 
pinning, the high-drag limit has been implicitly assumed through the following relationship between ft and /?': 

P = i-cL=Ypx?=np, (is) 

where 1Z is a dimensionless drag coefficient. In this drag description, imperfect pinning (a « 1, f3 « 1) corresponds to 
1Z » 1 so that eq. (|18|l requires /3 » a. Eq. (118[) is not true in general. The presence of non-dissipative forces between 
vortices and the solid to which they ar e pinned ca n give j3 « a for a and /3 both small, a regime of low drag that does not 



follow from eq. (|18p for any value of 1Z ( Link 2009) . As discussed below, it is the low-drag regime that is likely to be realized, 



with vortex creep being unstable in this regime. A crucial feature of our analysis is that we do not assume eq. (|18[) . 

To examine the stability of superfluid flow with imperfect pinning, we use a local plane wave analysis in the frame 
rotating with the normal matter at angular velocity tt n , in which v n — and the unperturbed flow velocity arising from spin 
down of the crust is Vo. Restricting the analysis to the regime A: AT? >> 1, where AT? is the thickness of the inner crust, the 
background flow can be taken to be uniform and the local analysis is valid. We take the unperturbed vortex lattice to be 
locally undeformed (cr e i = 0). The unperturbed creep velocity in the rotating frame (v„ — 0) follows from eq. (|17[) : 

— — - = av - puj X v - (19) 
at 

Below we estimate dr v o/dt ~ 10 -5 vo for a typical neutron star. Linearizing eqs. (|12|l and (|13p about vq, dr v o/dt, and er e i — 0, 
and neglecting dr v o/dt compared to vq, gives 

V ■ Sv = (20) 
^ + vq ■ VSv + 2U n x Sv = -VSp ~ X7S<j> - a el /p + Sf/p (21) 

p2fl„ x tSv J + p5u> x vo = -(Tei + Sf (22) 

where S denotes a perturbed quantity, and // = p, — p(Sl n x T ') 2 /2. We assume that a and /? are constants. The perturbed 
force is then 

S f = (1 - a)p5 {u: x v} + PpS {w x («xt))} + (l-7)(T e i(4 (23) 

The vorticity appearing in this equation is the total vorticity evaluated in the laboratory frame. For Vx»o << 2fl n , a good 
approximation for most neutron stars, the vorticity is 

u> = 2f2 n + V x Sv. (24) 

The final term in eq. (|23[) . associated with stress in the vortex lattice, will turn out to be negligible for vortex creep driven 
by a flow Vq » ct and vq » cv- 

We take the rotation axis to be z, with the unperturbed flow in the azimuthal direction, and along x at some point. For 
simplicity, we restrict fc to lie in the x — z plane, with an angle 6 with respect to the rotation axis. We further restrict the 
analysis to the quadrant 9 ^ 7i"/2. For shear perturbations, fc • Sv = 0, that is, the velocity perturbations in the directions 
y and e = — cos 9 x + sin 9 z are orthogonal to fc. 

We now Fourier transform (oc e lk ' r - l<Tt ) e q S (|20|) - (|23| ) and take the projections onto y and e. Defining a' = a — kvo s'm9, 
c = cos 9, and s = sin9, we obtain the system of equations: 
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—iff' + 2fl„f3 — i(l — a)kvos 


-2fin«c 7(ct s2 + 4c 2 )fc 2 







y 


Sv 


—iflkvosc + 2fi n ac 


-icr' + 2fi n /?c 2 - i(l - a)fcu s 


/ 2 2 2 2\ t 2 




e 


Sv 


— i(ff' + kvos) 





i2fi n o"'/c 




y 


r v 





— i(ff' + kvos) —i2£l n co' 







e 


r v 



The resulting dispersion relation is quadratic: 
(a 1 ) 4 + a3((j') 3 + a2(V) 2 + Oiff' + do = 0, 
where, in units with fi n — 1, 
a 3 = 2(1 - a)s{kv } + 2i(l + c 2 )/3 



o 2 = (a - l) 2 s 2 {kv } 2 + (2ips(l + c) - 2ia/3s + {kc T } 2 + ^■ysc 2 {kc v } 2 ^j {kv } 



4(a 2 + p 2 )<? — a r ys i {kcT} 2 — «7C 2 (1 + c 2 ){kcv} 2 + -'y 2 s 4 {kcT} i — -7 2 c 4 {fccv} 4 



= 0(25) 

(26) 
(27) 

(28) 



ai = ij/?7({fccT} 2 s 4 + {kcv} 2 c 2 s 2 ){kvo} 2 — (^aj[{kcT} 2 s 5 + {fccv} 2 sc 2 (l + c 2 )] — ^7 2 [{fccT} 4 s 5 — {kcv} 4 sc 4 ]^ {kvo} (29) 

a = jj 2 ({kc T } 4 s a - {kcv} 4 s 2 c 4 ){kv } 2 , (30) 
where, {fc«o} = fcvo/fin, {kcr} = kcr/^n, and {fccy} = fccy/fin- 



4 WAVE SOLUTIONS WITHOUT FLOW 



Before turning to the full problem with non-zero Vo, we consider the limit of vo = for the two cases of zero pinning and 
imperfect pinning. The dispersion relation is quadratic 

a 2 + 2ip{l + c)o - 4c 2 (a 2 + /3 2 ) + jc 2 T k 2 s 4 (-a + \lc 2 T k 2 ^ - -yc 2 v k 2 c 2 (a(l + c) + i-yc 2 v k 2 c 2 ^ = 0, (31) 
For zero pinning force (a = 7 = 1, /3 = 0), the dispersion relation to order c\k 2 and Cyk 2 becomes 



a 2 = (2fi n cos 8) 2 + c v (k 2 cos 2 9) (1 + cos 2 9) + c T k 2 sin 4 9, (32) 

as found bv lBavm fc Chandler ( 19831 ). The fluid supports Tkachenko modes for 9 — w/2, and axial modes (modified inertial 
modes) for 9 = 0. In the limit ct = cv = 0, the system supports only ordinary inertial modes. 

The role of pinning can be seen by considering axial modes for the case 7c 2 / k 2 « afi 2 . The solutions to eq. (I3ip in this 
limit are 

i74fc 2 ~ 



a± = 2i/3 fi„ ± 2a fi n + 



(33) 



which shows the damping effect of j3. Pinning strongly suppresses the axial mode given by eq. (|32[) . eliminating it entirely for 
perfect pinning. The waves are underdamped for /3 < a, which defines the regime of low drag that we will study further. 



5 INSTABILITY 

We now show that a non-zero background flow vq drives a hydrodynamic instability if the vortices are imperfectly pinned 
(a << 1, f3 << 1, 7 << 1). We are interested in flow velocities of vo ~ 10 5 cm s _1 . By comparison, 



c v ~ 10 c T ~ W~ v 



(34) 



We will find that there is an instability for wavenumbers k > fi n / vo . The hydrodynamic limit imposes the restriction kcr /fin < < 
1. The regime of interest is thus, 



fin/«0 < k « VL n /c T - 



(35) 



We will estimate below that a ~ 10 (3 ~ 10~ 10 . We assume that 7 is similarly small. 

We will not present here an analysis of the full mode structure of the system, but focus on two low-frequency modes 
that appear for imperfect pinning. We simplify the problem by proceeding to linear order in the small quantities kcr /fin and 
fccy/fi„. At this level of approximation: 



a 3 = 2(1 - a)s{kv } + 2i(l + c)/3 



(36) 
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a 2 = (a - l) 2 s 2 {kv } 2 + 2i/3s{l + c 2 - a){kv } - Ac 2 {a + (3 2 ) (37) 

ai = «o = 0, (38) 

that is, the vortex lattice exerts no stresses on the fluid to first order in ct/vq and cy/vo- The dispersion relation eq. (1261) 
now simplifies to: 

(a') 2 {(cr') 2 + 2({1 - a}kv s + iP{l + c 2 }) a' + (1 - a) 2 k 2 v 2 s 2 - Ac 2 (a 2 + /3 2 ) + 2i/3kv s(l + c 2 - a)} = 0. (39) 

The two degenerate solutions (cr') 2 = correspond to a — kvo sind = k ■ vg, the frequency associated with translation of the 
wave pattern at velocity vq. The other two other solutions are, switching from a' to a and restoring Q n , 

a± = akvo sind — iQ n (l + cos 2 8)f3 ± (4f2 2 a 2 cos 2 8 — fi 2 /3 2 sin 4 8 — 2ia/3Q„kvo cos 2 #sin#) 1//2 . (40) 

For /3 << a and low wavenumber kvo « &n, there are two damped modes 

o± ~ ct(kvo sin 8 ± 2f2 n cos 8) — i/3 ^0, n {l + cos 2 9} ± t^ u o cos 8 sin . (41) 

Slow vortex motion has introduced two low- frequency modes to the system. Removing pinning and drag (a = 1, /3 = 0) and 
taking k = 0, we recover the ordinary inertial modes a± — ±2Q, n cos 8. 

Above a critical wavenumber k c , the the solution with eigenvalue er_ is unstable: 

k>k^2^ 2 + a ^ /2l + COs2e . (42) 
vo a sin 8 cos 6 1 



Numerical solution of the full dispersion relation, eq. (I25p . for reasonable values of ct, cy, and vo, confirms that there are 
no other instabilities. The critical wavenumber k c is minimized for 8 — tan _1 (\/2). For k » k c , we have the approximate 
solutions 

a± ~ akvo sind =p i(a/3 Q.„kvo cos 2 8 smd) 1 ^ 2 . (43) 



The instability arises from coupling between velocity and vorticity through the first two terms of eq. (|23[) . Dissipation 
damps perturbations for k < k c , but for k > k c the finite vortex mobility gives rise to growing perturbations under the Magnus 
force. For k >> k c , the growth rate scales as (aPvo) 1 ^ 2 - For << a, k c takes a constant value, but the growth rate of the 
mode becomes small, going to zero as /3 goes to zero. In the highly-damped regime, ft » a, damping restricts the unstable 
mode to large k, generally stabilizing the system. There are no unstable modes for either a = or ft = 0; the instability 
occurs only if the vortices move with respect to the crust, both along the flow and transverse to the flow. 

We now show that our neglect of shear deformations of the crust is a good approximation. The modes we are studying 
are in the regime kvo > fin- In this limit, the dominant contribution to the shear force per unit volume in the fluid is (see eq. 

EU) 

8f/p ~ kv 8v. (44) 

Because the vortices are nearly perfectly pinned, this shear force creates a strain field in the solid with a shear force per unit 
volume of 

Sfs/p ~ c 2 k 2 Su, (45) 

where c s is the shear speed of the solid and 8u is the characteristic displacement. The speed of a mass element in the solid is 
Sv n ~ Ke(a±)Su. Equating Sf and 5f s , and using eq. (I43|) for the limit of large k, gives 

Sv 

The values a = 10" 10 , v = 10 5 cm s 1 , and c s = 10 s cm s , give 5v n /5v ~ 10 16 . The displacement of the solid is very 
small for two reasons: i) the solid is very rigid compared to the vortex lattice, and, ii) the vortex creep modes are of very low 
frequency, proportional to a « 1. 

We have restricted the analysis to shear waves (V • 8v = 0). Because these waves to not perturb the den sity, we do no t 



expect the finite compressibility of the matter to change our results. Compressibility will introduce new modes (Haskell 2011) 
an effect that merits further study in the context of imperfect pinning. 



6 ESTIMATES 

To obtain the growth rate of the instability, we now estimate the pinning parameters a and /3 for the vortex creep process. To 
make these estimates, we regard the process of vortex creep as consisting of two distinct states of motion for a given vortex 
segment. Most of the time, the vortex segment is pinned. A small fraction of the time, the vortex segment is translating 
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against a drag force to a new pinning configuration. The mutual friction force we are using feq. U6[) is, ignoring the small force 
from the vortex lattice, 



f/p 



- aui x v + /3tl> x (w X v). 



(47) 



This force represents the average force exerted on the neutron fluid by the vortex array. The first term is the Magnus force for 
perfect pinning, while the remaining terms give the contribution to the force due to vortex motion. For those vortex segments 
that are unpinned and moving against drag, we take the force to have the same form, but with different coefficients: 



fo/P = 



ao x v + /3o uj x (w x v). 



An unpinned vortex segment remains unpinned for a time to ~ d/vo, 



(48) 

where d is the distanc e the segment moves before 
~ 19931), 



repinning. This distance is comparable to the distance between pinning sites (jLink et al.lll993h . roughly ten times the unit 
cell size, giving to ~ 10 -15 s, much shorter than the hydrodynamic timescales of interest. Suppose that at any instant, the 
fraction of vortex length that is unpinned is f v « 1. We now average f over a volume that contains many vortices, and 
over a time long compared to to but short compared to hydrodynamic timescales, to obtain 



(fo/p) — w x v ~ fv a o w x v + fvPa w x (w x v). 



(49) 



Quantities related to the flow are unchanged by the averaging procedure since the superfluid flow velocity is independent of 
whether vortices are pinned or not. The factors of f v in eg. (149 II account for the fact that only the motion of the translating 
vortex segments contributes to the mutual friction (see, also, Jahan-M iri 2006). The value of /„ is unimportant for the following 
estimates. 

The force of eq. (|47[1 . which is appropriate for vortex creep, must equal the average force {fo/p), giving the following 
relationships: 

a = f v a and = f v #> => £ = ^ (50) 

a ao 

We now use estimates of fio/ao to obtain the ratio /3/a. 

The dominant drag process on unpinned vortex segments considered so far arises from the excitation of Kelvin modes 
as the vortex moves past nuclei. Calculations of dissipation by Kelvin phonon productio n on a long vortex with periodic 
boundary conditions for vo ~ 10 r cm s" 1 give typical values of Po/o-o = 0.1 and qq ~ 1 ( Epstein fe Bavm 1992). Pinning 
occurs for «oi; 10 5 cm s _1 , a nd Po/ap is likely to be significantly smaller in this velocity regime due to strong suppression 
of Kelvin phonon production ([JonesJ|l992j). Vortex creep is therefore a low-drag process if Kelvin phonon production is the 
dominant dissipative mechanism. We fix j3/a = 0.1 for illustration in the following, which we consider to be an upper limit; 
we expect typical values to be smaller. 

We now estimate ft. We adopt polar coordinates (r,<j>,z), with the unperturbed vorticity along z and the unperturbed 
flow Do along </>, and take the unperturbed flow and vortex velocity field to be axisymmetric. In the rotating frame, the 
unperturbed vortex velocity from eq. (|19[1 is 



dr v0 



= avo<j) + Pvof = h 



dr v0 



(51) 



dt ^ ^ dt 

where n is the average direction of vortex motion. 

For steady spin down of the star, the inner crust superfluid and the crust are spinning dow n at the same rate for a local 



differential velocity Vq. The creep velocity in this steady state is related to the spin-down rate by (|Alpar et al.ll 19841 : iLink et al 
19931 ) 



Q = -2 



dr v0 
r dt 



r = —2— vo p = fio, 
r 



(52) 



where SI is the spin rate of the superfluid, Slo is the observed spin down rate of the crust, and r is approximately the stellar 
radius R. We arrive at the estimate 

-i 



P = 



R 



4v t a 



10 



"( 



vo 



10 5 cm s 



( U 



\ 10 4 yr 



(53) 



where Q ~ Qo is assumed, and £ age = r2o/2| J7o I is the spin-down age. Eq. (|53|l . with /3 — 0.1a, gives the fiducial value 
a/3 = 10~ 21 . For this value, we deduce /„ ~ (a/3/aoPo) 1 ^ 2 ~ 10 -11 , that is, most of the vortex length is pinned at any 
instant. The unperturbed vortex creep speed, from eq. (|51[) . is ~ avo ~ 10 -5 cm s _1 << vo, justifying the neglect of dr v o/dt 
compared to vo in the stability analysis. 

We can now proceed with estimates of the instability length scale and growth rate. For /3 < a and 8 — tan - (v2) m eq. 
()42|) . the critical wavenumber is 

n \ / wo 



6— = 6 x 10~ 3 

vo 



( " ) 

V100 rad s -1 / 



10 5 



(54) 
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corresponding to a wavelength A = 2-n /k ~ 10 m. For k >> k c , the growth rate from eq. (|43|l is 

_ Mff _ )=i 0.6^J i 100 rads-0 U cm s~0 It^J yr ' (55) 
The hydrodynamic treatment is restricted to kcr << Q. To estimate how high the growth rate could be, we consider a 
maximum wavenumber defined by CTfc max = 0.1 Q, where ct — 10 _1 (£7/100 rad s -1 ) 1 / 2 cm s _1 . The growth rate at this 
wavenumber, from eq. ()43[) . is 



For = 100 rad s~\ the corresponding wavenumber is fc max ~ 100 cm -1 . Eq. ()56|) does not represent a physical limit, but 
only the restrictions of the hydrodynamic treatment; the instability could continue to exist also for wavenumbers in the regime 
kc T > n. 

If vortex creep is in the strongly-damped regime /3 » a, contrary to the estimates here, there is still a broad window 
for instability. Requiring k c < k max fflVGS 

P < 2 x 10 4 ( — T ) ( ) ^ a , (57) 

V10 5 cm s -1 / V100 rad s" 1 / 

and the star will be unstable at some wavenumber that is consistent with the hydrodynamic regime kcr << O. 



7 DISCUSSION AND CONCLUSIONS 



We have identified a dissipation-driven instability that could operate in the fluid of the neutron star crust over length scales 
shorter than ~ 10 m and over timescales as fast as months. This instability is different than other superfluid instabilities 
already considered in two-component systems. In the case of liquid helium, the Glaberson-Do nnelly instability arise s when 



the fluid normal component has a component of flow along the rotation axis of the system (jGlaberson et al.lll974f h That 



instability occurs even if the mutual friction is zero. A variant of the Glaberson-Donnelly instability has been studied in 
the mixture of superfluid neutrons and superconducting protons of the neutron star core, and instability was found to occur 
under the assumption that the vortices are perfect l y pin n ed to the flux tubes that p e netra te the superconducting proton 
fluid ( Glampedakis et al. 20081 : van Hoven fc Levin 20081 ) . Glampedakis fc Anderssonl (|20Qgh have identified a two-stream 
instability that might occur in the neutron-proton mixture of the core, again assuming perfect pinning of vortices to flux 
tubes and neglecting magnetic stresses. By contrast, the instability described here exists in a single-component fluid, and 
occurs because the vortices can move with respect to the solid. The component of the motion that is transverse to the flow, 
the component related to /3 in eq. (|16p . though dissipative, is essential for the instability to occur. 

To illustrate the basic instability, we have taken the pinning coefficients a and j3 to be constants. For thermally-activated 
vortex creep, these coefficients will have exponential dependence on the velocity difference between the superfluid and the crust 
( Alpar et al.ill984l ; iLink et al.|[l993T ) . We expect that this strong velocity dependence will significantly enhance the growth rate 
of the instability. The instability of the system will be determined by four coefficients: a(vo), /3(vo), and the derivatives da/dv 
and d/3/dv evaluated at Vq. Further work is needed to calculate these coefficients and to incorporate them in the stability 
analysis. We restricted the analysis to wavevectors that are co-planar with the rotation axis and the unperturbed flow; more 
general perturbations should be studied, including waves with compressive components. 

The flow of the inner crust superfluid could become turbulent. If the system evolves into a state of full y-developed super- 
fluid turbulence wit h a vortex tangle, the f riction force in the fluid would be better described by an isotropic (jGorter fc Mellink 
1949) or polarized l|Andersson et al.l l2007) form, rather than the anisotropic form of eq. (|16|l that is appropriate to a regular 
vortex array in the initial stages of the instability. The friction force for the tangle could be larger or smaller than the force 
given by eq. (|16|l . On the one hand, the tangle has more vortex length per unit volume to interact with the solid, which tends 
to increase the force for a given value of the velocity. On the other hand, if the vortex distribution becomes highly tangled 
the momentum transfer from different regions will cancel to at least some extent, decreasing the friction force compared to 
that of a straight vortex array. If the force increases, so does the effective value of /?, and the average value of the equilibrium 
differential velocity will decrease (see eq. 153 p . This effect could spell trouble for inner-crust models of glitches. By contrast, 
if the friction is decreased by turbulence, there will be more excess angular momentum in the superfluid available to drive 
glitches. 

The instability results from the forcing of the vortex lattice through the solid lattice. If fully-developed turbulence results , 
the closest experimental analogue might be grid turbulence that has been well studied in superfluid helium ( Smith et al.lll993l ). 
Further work is needed to determine the observational consequences of this instability, or if it is damped by some mechanism 
that has been overlooked here. We have sh own elsewhere that a similar instability occurs in the mixture of superfluid neutrons 
and protons of the outer core |Linkll2012l ). 
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